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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 

ADVANCE RESTRICTED REPORT - 

COMPARISONS OF METHODS OF COMPUTING BENDING 
MOMENTS IN HELICOPTER ROTOR BLADES IN THE 
PLANE OF FLAPPING. 

By John E. Du berg and Arthur R. Luecker 
SUMMARY 


Several methods of computing bending moments in 
helicopter rotor blades In the plane of flapping are 
reviewed, and the results of a numerical example analyzed 
by four different methods are compared. The effect on 
computed bending moments of the tip-loss correction 
introduced by Wheatley is also considered. 

The comparison indicates that, from the standpoint 
of accuracy of results and esse of application, the 
method proposed by Cierva is the most suitable for routine 
analyses. The tip-loss correction Is shown to have a 
substantial effect on computed bending moments. 


INTRODUCTION 


The determination of the bending stresses In rotor 
blades during flight Is one of the important problems In 
the stress analysis of the structure of the helicopter. 

The problem is complicated by the fact that the air loads 
and the Inertia loads are continuously changing as the 
blades rotate and that the bending deflection of the 
blades has an important effect on the moments developed 
therein. The analysis can be simplified if the blades 
are assumed rigid and therefore unable to bend, but under 
this' assumption the computed bending stresses superimposed 
on the centrifugal-tension stresses are relatively high. 

A more exact analysis, in which the relieving effects of 
centrifugal tension on blade bending are included, gives 
much lower values for these calculated bending stresses. 

In the present paper comparisons are made between 
the various methods that h?ve been proposed for the 
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calculation of the bending stresses in the plane of 
flopping. The effect of introducing a correction for 
tip loss # which has heretofore been ignored in stress 
analyses, is also considered. A numerical example is 
analyzed by each of the several methods to provide a 
concrete basis for the comparisons. 


THEORY AND BASIC ASSUMPTIONS OP METHODS OP ANALYSIS 


The basic rotor theory on which stress analyses of 
rotor blades are based is that of Glauert, Lock, and 
’.Vheatley. (See references 1, 2, and 3# respectively.) 

A consideration of the forces acting on the rotor blade 
during steady forward flight, as given by this theory, 
indicates that the total bending moment at any radial 
station may be resolved into the following three components 

(1) A component that is identical with the bending 

moment in a rigid blade and is therefore 
Independent of blade bending 

(2) A component, due to axial tension, that depends 

on the blade deflection and therefore varies 
as the blade bends 

(3) A component due to the inertia forces associated 

with the variation of the deflection with time 

In the simplified theories. In viilch the blade is assumed 
static and rigid, the second and third components do not 
occur. T -?hen flexibility Is taken Into account so that 
the second component Is Included, the problem Is compli- 
cated because of the Interaction between the moments and 
the blade deflections. If the third component Is included 
the problem Is further complicated because the elastic 
curve must satisfy dynamic as well as static conditions. 


METHODS OF ANALYSIS 


The methods of analysis for rotor-blade bending 
moments in tho- plane of flapping that are available In 
the literature fall into three categories as follows: 

(1) "Exact 1 11 analyses In which Inertia forces due to 

rate of change of blade deflections, as well 
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as the Interaction between blade deflections 
and moments, are considered 

(2) Analyses, based onstatic loadings, that Include 

Interaction between blade bending and moments 
but neglect Inertia forces due to blade- 
bending deflections 

(3) Simplified analyses In which the bending moments 

for the rigid biade are computed and approxi- 
mate corrections to account for the reduction 
In the bending moment due to the axial load 
are applied; the reduction due to the axial 
load has sometimes been called centrifugal 
relief 

The analysis given In reference Ij. Is believed to be 
the only published exact analysis and Is based on the 
rotor theory of references 1, 2, and 5» This theory has 
been refined and extended In roference 3« In appendix A 
of the present peoer, one of these refinements - allowance 
for the reduction In lift near the blade tip - Is Intro- 
duced Into the analysis of blade bending. 

Most of the solutions that have been presented In the 
literature fall Into the second category, in which the 
problem treated Is purely a static one. This approach 
reduces the problem to that of a beam under combined 
bending and tension and, therefore, the methods developed 
for such problems can be applied to the analysis of rotor- 
blade moments. In reference 6, a solution of this class, 
which makes use of "type" solutions to facilitate the 
numerical calculations, is given as a simplification of 
the analysis of reference Ij.. Two other methods of solving 
the static problem are given in references 7 and 8. In 
reference 7 the solution is effected by means of the theo- 
rem of three moments generalized to include the effect of 
centrifugal tension. In reforenco 8 the solution Is 
obtained by finding a deflected shape for the blade that 
is consistent with a minimum of potential energy for the 
system. 

The methods of the third category, which should have 
more appeal for use in practical applications, provide 
simplified formulas for obtaining the bending moments In 
the blade by correcting the rigid-blade bending moments 
for the relief caused by centrifugal tension. Cierva has 
proposed a formula, given in reference 6, which states 
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that the moment developed In the actual rotor hlade is 
equal to the product of the rigid-blade bending moment 
and the perfectly flexible blade bending moment divided 
by the sum of these two moments. The perfectly flexible 
blade bending moment Is a fictitious bending moment 
obtained by multiplying the stiffness of the actual blade 
by the curvature of the perfectly flexible blade. The 
solutions for the rigid-blade bending moment and the 
flexible -blade bending moment are given in appendix B. 
Hohenemser in reference 9 has developed a method of com- 
puting the blade bending moment by multiplying the rigid- 
blade moment by a correction factor, given herein in 
appendix B, that depends on the blade radius, blade 
stiffness, and centrifugal force at the root. The formula 
is developed on the basis of uniform distribution along 
the length of the blade of both the blade mass and blade 
bending stiffness. Reference 10 contains a brief discus- 
sion of Hohenemser* s formula as well as the development 
of formulas for the calculation of the rigid-blade bending 
moments. 

Although wind-tunnel and fligit tests have proved 
the basic rotor theory reliable when applied to the rotor 
as a whole, the theory contains approximations that raise 
some doubts as to the order of accuracy in computing air 
forces at particular points in the rotor disk and, hence, 
doubts as to the accuracy of bending moments computed by 
any of the methods mentioned herein. These doubts can be 
removed only by further tests in which blade bending is 
measured under operating conditions. 


COMPARISON OP BENDING MOMENTS COMPUTED BY 
DIFFERENT METHODS 


In order to compare the bending moments obtained by 
the several methods of analysis , calculations were made 
for a rotor with three blades, each having the following 
physical properties: 


Radius, feet ... 

Mass, slug per foot ....... 

Blade chord, inches 

Bending stiffness, El, pound-feet^ 
Pitch setting (untwisted), degrees 


12.5 

O.0519 

• 4 

76^0 
. 10 
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The rotor was assumed to be rotating at 370 rpm and to 
have a forward velocity of 100 miles per hour, which 
corresponds to a tip-speed ratio p. of 0.30. The ratio \ 
of axial velocity to tip speed was assumed to be -0.079? 
the negative sign indicates flow downward through the 
rotor disk. 

The bending moments In the plane of flapping have 
been computed for the rotor blade described in the 
preceding paragraph by methods in each of the three 
oatagorles of analysis mentioned in the preceding section. 
The moments computed for four azimuth angles (i|/) are pre- 
sented in figure 1. Inspection of these moment diagrams 
reveals that only a very small difference in computed 
bending moment at any station exists between the values 
given by the exact method and this exact method modified 
by neglecting the inertia forces due to blade bending. 

The small magnitude of this difference substantiates the 
assumption, common to all approximate methods of analysis, 
that these Inertia forces are negligible. 

The method of analysis proposed by Gierva, which Is 
given in detail in appendix B, gives a maximum bending 
moment for any azimuth angle that differs only slightly 
from the maximum bending moment given by the exact method 
for the same azimuth angle. The moment diagrams given by 
the Cierva method are In close agreement with the moment 
diagrams for the exact method except for the outer third 
of the blade. In the outer third of the blade the bending 
moments predicted by the Cierva formula decrease to zero 
less rapidly than those predicted by the exact method, but 
this fact Is of little consequence because In this part 
of the blade the moments are less than the maximum moment. 

The method of analysis proposed by Hohenemser gives, 
as shown In figure 1, a bending moment diagram at each 
azimuth angle that differs appreciably from those given 
by the other methods. The computed moment Is higher In 
the inner part of the blade and lower in the outer part 
as compared with the moment computed by the exact method. 
At an azimuth angle \| / of $ 0 °, where the maximum bending 
moment approaches its smallest value, the Hohenemser 
method gives a maximum bending moment appreciably higher 
than that given by the exact method. At \|r - 270°, where 

the maximum bending moment approaches its largest value, 
the Hohenemser method gives a maximum bending moment 
appreciably smaller than that given by the exact method. 
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The labor required to evaluate the bending stress 
by the approximate methods of Clerva and Hohenemser Is 
only a small fraction of the labor required by the more 
exact methods. Of these two approximate methods, that of 
Clerva is recommended for practical stress analysis of 
helicopter blades because the bending stresses obtained 
thereby are In good agreement with those given by the exact 
method. 

In figure 2 are shown the results obtained by the 
exact analysis extended to include the effect of tip loss 
by the method introduced by Wheatley (reference 3)« The 
tip-loss factor was assumed to be 0 . 97 » which means that 
the outer 3 percent of the blade is assumed to produce 
no lift. The method of analysis used assumes zero moment 
at the 0.97 point and does not consider bending between 
this point and the tip. At all azimuth angles, the 
consideration of tip loss reduces the maximum bending 
moment considerably and causes a reversal of bending near 
the tip. Because the assumption made concerning tip 
effect is only a crude approximation of the actual lift 
distribution in this region, the calculated moments near 
the tip may be very much In error. The large reduction 
in the bending moment near the center, however, Is 
significant. 

The bending stresses at each station In the rotor 
blades of a helicopter vary between a maximum and a 
minimum with each revolution of the rotor. For the par- 
ticular blade analyzed, which was assumed to be of all- 
metal aluminum- alloy construction with a section modulus 
of O .167 Inch cubed and an effective cross-sectional area 
of O .56 square Inch, the total stress, which includes the 
centrifugal tension and the superposed bending stress, 
varied between the limits shown In figure 3 for the upper 
and lower fibers at the different stations along the 
blade. Figure 3 is based on the results of the exact 
bending-moment analysis, neglecting tip loss. 


CONCLUSIONS 


A comparison was made of rotor blade bending-moments 
obtained by several methods of analysis. The results for 
the single numerical exairple studied Indicate that the 
additional accuracy obtainable by the "exact 11 analysis, when 
compared with the best of the approximate methods now In 
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use, does not justify the extra labor involved in applying 
the exact method. 

Of the two approximate methods studied, that of Cierva 
gives the best over-all agreement with the more exact 
analysis and is well suited for the determination of 
stresses in blades. The method proposed by Hohenemser, 
although easier to apply than the Cierva method, is less 
desirable because it gives maximum bending moments that 
may differ by as much as one-quarter from those computed 
by the exact method. 

Consideration of tip loss results in a substantial 
reduction in maximum bending moments and should receive 
further attention in future studies of methods of blade 
stress analysis. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. , Jui 
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APPENDIX A 

EXACT METHOD OP CALCULATING ROTOR -HLADE BENDING MOMENTS 
Differential Equation of Blade Bending 


The basic assumptions of the blade-bending analysis 
referred to herein as the exact method are identical with 
those used in reference I4. for the development of the dif- 
ferential equation except that a correction is added to 
include the effect of tip loss. The correction for tip 
loss follows the method developed by Wheatley in refer- 
ence 3 in which the air load on the outer few percent of 
the blade is neglected. 


Figure Ij. shows the coordinate system used and the 
nomenclature involved in defining the deviation of the 
flexible blade from the rigid blade. The symbols used in 
the analysis are defined in appendix C. In figure 5 a **e 
shown the forces and moments acting on each blade element. 
A consideration of the equilibrium of the element in the 
direction parallel to the longitudinal axis of the blade 
and in the direction normal to the blade yields the 
following differential equation for the blade bending 
deflections : 



(> 


- K [1 + 


_2QRq 

mB 2 R 2 




dx e 


ZKX sz + ss iz. _ iisit 

■ . n 2 dt 2 EI dr 


(Al) 


If the blade mass Is uniformly distributed from hinge to 
tip, the mass of the tip section, which is assumed to be 
without air load. Is 


Q = mR(l - B) 

and the distance from the tip mass to the hinge Is 


Rq = |(1 + B) 


If these relationships are substituted in the second term 
of the differential equation some simplification results 
and the differential equation becomes 



* 2 \S5Z + 2k £Z + 22 5S 

/ dx 2 dx n 2 dt 2 EI dr 
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If the tip loss is neglected, B reduces to unity 
and the -differential equation reduces to that given In 
reference I4.. 


Blade Loading 

The term on the right-hand side of equation (Al) 
represents the blade load gradient, which is assumed 
to be the same as that on a rigid blade. The blade 
loading can be evaluated by means of the rotor theory 
given in references 1 and 3- In the present analysis the 
blades are assumed to be built with constant chord and 
without twist and the effect of periodic twist of the 
blades is neglocted. The effect of both types of twist 
can be considered by the methods of references 2 and 3» 

If the instantaneous flapping angle of the blades Is 
represented by the equation 

p = a Q - ax cos - b^ sin \|/ ( A2) 


the 

the 


The 


velocity components causing aerodynamic forces on 
blades are 

Ut = xBRfl + fiRfl sin \|/ 

Up = \Rfi - xERft(a;L sin vU - b^ cos U') 

-nRfl(a 0 - a^ cos \|/ - b^ sin v)cos \| / 


>U3) 


load gradient Is then given by 




dS 

dr 


— “cpa(uipUp + 0Uip2^ - mBR^a 0 x - mg 


(Al+) 


Introduction of the relationships given in equa>- 
tlone (A3) into equation (AJ+), expanding, and neglecting 
higher harmonic e gives for the load gradient 
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dS 

dr 


|opan^n 2 sb 2 * 2 + b(\ - + |ep 2 

+ £ -a^B 2 * 2 + 20 jj,Bx + p,\ + s ^ n ^ 

+ £b]B 2 x 2 “ a oM-B^ + ^ji 2 bi^J co 3 \J/^ 


jpaR 2 ^ 2 


(A5) 


Flapping Coefficients 

The flapping coefficients a Q , a j_, and b 1# which 
define the blade motion given in equation (A2) , are deter- 
mined from the condition that the moment at the flapping 
hinge is zero for all azimuth angles. The coefficients 
are : 



where 



i! = QERRq + ^mB3R3 
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Solution of the Differential Equation of Blade Bending 

The method of solution ' of ' the ‘differential- equation, , 
follows that given in reference Ij., The blade deflections 
are assumed to be expressible in the form 

7 = 7i + 72 sln ^ + cos \Jf (A7) 

in which y-^, * ancl are functions of the distance 

along the blade span. Substitution of equation (A7) and 
(A5) in the differential equation (A3) yields th9 following 
differential equations for the functions yi, and 7$* 




The four boundary conditions that the functions yi , yg, 
and y^ must satisfy are : At the hinge the deflection 
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and moment are zero, and at radius HR the moment is zero 
and the shpar is that required to maintain equilibrium of 
the tip mass Q. The last-mentioned boundary condition, 
at x = 1.0, can be expressed in the following form: 

dy l _ 

— + BRa 0 + 
cx 




d'72 dy 2 
A *^J- = Hx- ‘ ^ 


A— = 
dx^ 


d 7 7 

ax - 7 3 


where 


El 

A = 

QD^rSRo 


The differential equations for yi, y 2 , and y^ can 
be solved approximately by the method of " collocation" . 
(See reference li.. } The method of collocation consists 
essentially of expressing the solution as a linear com- 
bination of functions that satisfies the boundary condi- 
tions indeoendently of the value of the coefficient 
associated with each function. The combination of 
functions is substituted into the differential equation, 
and the ccefficienta are so adjusted that the resultant 
expression satisfies the differential equation at as 
many points as there are functions. 

When the tip effect is considered the following 
set of functions can be used: 



1 
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y 3 



L_±_1 x p+3N 
P + 3 J 


Lli x p+3) 
P + 3 J 


If the tip effect is neglected, the form of the 
solution for y^ Is 


71 = 


°ll x 


n 

\ 

/_ 

P=2 



_!L_ x p+2 

p + 2 


+ P(P + 1? 

(p + 2 ) (p + 3 ) 


The form of solution for y 2 and y* remains as 
before. 5 


Numerical Example 

A numerical example Is presented with tip loss 
neglected, for a blade with the following properties 
and operating conditions: 

R = 12.5 feet 

c = 9~ inches 
2 

m = 0.0519 slug per foot 
a = 5.73 

El = 76 I 1 .O pound-feet 2 

= 38.8 radians per second 
ix = 0.300 
\ = -0,079 

6 = 10 ° = 0.175 radian 
p = 0.00230 slug per cubic foot 



NACA ARR No. L5E23 


llj. 

If the tip loss is neglected, 

B = 1.0 and Q = 0 

Equations (a 6) give the following values for the 
flapping coefficients: 

a Q = 0.077922 

a x = O.O96963 
= 0.029827 

The factors y and K have the following values: 

y = 7.5535 
K = 121^.83 

If six functions are chosen for the solution to the 
differential equation for yj, the form of the equation is 

n = on* + =i2 (* 3 - ^ + 15* 5 ) + c u - 1* 5 + 1* 6 ) 

+ M ♦ |2x?) ♦ o 15 - ft + 15 x 8 ) 

+ 0 l 6(* 7 - |* 8 + (AH) 

When this function Is substituted in the differential 
equation (a 8), the equation Is satisfied at the points 
x = 0, 0.2, O.lj., 0.6, 0.8, and 1.0 if 

Cu = -0,171 foot C]2^ = -0.919 foot 

C 12 = 0.133 foot = 0.569 foot 

= 1.197 feet = 0.503 foot 
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Substitution of the coefficients in equation (All) 
yields 

y 1 = -0.171X + 0.133x3 + 1.o61+x4- - 2.315x5 + 2.273x 6 
-0 . 74-8x7 - 0 .l 4 l 4 . 9 x 8 + 0.293x9 


The corresponding bending moments are obtained by 


FI 

differentiating twice and multiplying by =i f 

R 2 


which gives 


Mi = 14.8.9(0 .795* + 12.77X 2 - 14.6.30x3 + 68.18x4- 
-31. 4-0x5 - 25.17x6 + 21.12x7) 

By a similar process the results for M 2 and Mj corre- 
sponding to the deflections yg and yj are 

m 2 = 4-£.9(5.i6x - 142.19x2 + 136. 1x3 - 239.9x4- 

+ 232.3x5 - 109. 2x 6 + 17.8x7) 

Mj = 4-8.9(o .333* + 1.32x2 - 5.08x3 + 5.93*4- 

+ 1.59x5 - 9.18X 6 + 5.09x7) 


The general expression for the bending moments is 
M = Mi + M? sin \|/ + Mj cos \j/ 

At a distance from the flapping hinge of x = 0.6 the 
bending moment is 

M = 4-3*2 - 15*9 sin \|/ + 9*0 cos \|/ 
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The minimum bending moment at the station x = 0.6 
occurs when \|/ = 120° and is 

M = I4.3 *2 - 13.8 - J^.5 = 2 i |..9 pound-feet 

The maximum bending moment at the station x = 0.6 
occurs when \| / = 300° and is 

M = 1 +- 3 -2 + 13.8 + I4. . 5 = 61.5 pound-feet 
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APPENDIX B 

APPROXIMATE METHODS OP CALCULATING ROTOR- 
BLADE BENDING MOMENTS 
Cierva Method 


The Cierva formula for the bending moment in a 
rotor blade, given in reference 6, is 


M g Mf 
M r + 


(Bl) 


where M r is the moment in a rigid blade and Mf is a 
fictitious moment obtained by multiplying the curvature 
of a perfectly flexible blade by the actual blade bending 
stiffness El. 


If the blade Is assumed rigid and tip loss Is 
neglected, the rigid-blade bending moment obtained by 
Integrating the moment of the load gradient given In 
equation (A5) la 


M r 


■ •( - fa* - J 

[" 


+ j^C* + 20nB x + {mA + —^JAx 


sin 


f 


+ l b l c x “ M- a o B x + 


i^bi 

-jj— A x cos 'l' 


(B2) 


in which 


A x “ 




= l l tadx = K 1 - 21 + 

** = £ £ x6x6x = K 2 ■ 5x + * 5 ) 

G * = £ £ = ii ( 5 - ^ + 


> (B3) 




■ ■Ill II 


II II 


■ llll III 
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If the blade loading for the perfectly flexible 
blade la assumed equal to that of the rigid blade the 
fictitious moment existing in the blade Is 


M< 


= Sir 

9R 


f 


0 + - 


e - 


l8g 

YRfl 2 


-a x 


(1 + x ) 2 

+ ^aiu 2 


1 [i + 3n\ 


-a 7 - 


(1 + x ) 2 

b l + f b lM- 2 


sin \j/ 


h ■ 


(1 + xV 


cos \|r 


(bU) 


The numerical values of Mf and M r can then be 
combined to give the actual blade bending moment according 
to equation (Bl) . 

In reference 9 Hohenemser developed the following 
formula for the moment In a blade that has a uniform 
distribution of both blade mass and blade bending 
stiffness : 


M = 


Mr 


1 + 0.052- 


R 2 Pq 
El 


(B5) 


in which M r is the rigid-blade bending moment and p G 
Is the centrifugal tension at the flapping hinge. 

For a blade of uniform mass distribution, 

P Q = 

Substitution of the value for P c in equation (B5) 
reduces this equation to 


M = 


Mr 

1 + 0.052R 
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in which 

- - - _ mR^O 2 

2EI 


Numerical Examples 

The bending moments in the same blade and for the 
same operating conditions as were considered in appendix A 
are calculated by both Cierva’s and Hohenemser's method 
for x = 0.6 and \j/ = 120° and 300°. 


When the basic data from appendix A are substituted 
in equation (B2) the following result is obtained: 

= 192,000 £(0.1750C x - 0.1410B* + 0.00651A x ) 

- (0.0970C X - 0.1050^ + 0.0215A X ) sin t 

+ (o.0298c x - 0.023IJ.3X + 0.000 67 A x ) cos \|/J (b6) 

For x = 0.6, equations (B3) give 

A x = 0.0800 
Bx = 0.0693 
c x = 0.0608 

Substitution of- these values of A x , Bx, and C x in 
equation (B6) gives the following expression for the 
rigid-blade moment at x = 0.6: 


Jfp = 265 - 65 sin \j/ + I4.8 cos \|/ 


The perfectly flexible blade bending moment is 
obtained from equation ( Blj.) and is 


M f = 512 ([o.l750 - °' 19lt5 J + [-0.0970 + 0,1619 1 

II <1 + *> 2 J L (i + *> 2 J 


sin \|r 


+ 


( 


0.0298 


0.0318 
(1 + x) 2 
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At x = 0 . 6 , 

Mf = 50.7 - 17*4 sin \|/ + 8.9 cos 

At = 120 ° (the approximate position for minimum bending 
moment) , 

Mp = 265 - 56 - 2i|. = I85 pound-feet 

Mf = 50*7 " 15*1 “ Ij-.lj- = 31-2 pound-feet 

Therefore, the Cierva method (see equation ( B 1 ) ) 
gives for the bending moment M at the azimuth position 
\j/ = 120°, 

M = * 21 ‘Jr - 26.7 pound-feet 

105 +31 

Similarly, the Cierva method gives for ^ = 300 ° 

M = 59 . 14. pound-feet 

When Hohenemser ' s method is used, at \j/ = 120 °, 

M = 185 x = 21^.7 pound-feet 

1 + (0.052 x 12)+. 8) 

and similarly, at \J/ = 300°, 

M = 3 k -5 x = 14.6.0 pound-feet 


1 + (0.052 x 1214-.8) 
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APPENDIX C 


- - SYMBOLS 


A dimensionless coefficient 

a slope of lift curve 

a G constant term in Fourier series that 

expresses p 

a^ coefficient of cos i{r in expression for p 

bi coefficient of sin \|; In expression for p 

B tip-loss factor (blade elements outboard of 

radius BR are assumed to have no lift) 

c lpf C2p» Cjp coefficients in equations for y^, y£, and y^ 

c blade chord (constant) 

El flexural stiffness of blade 



e 


acceleration due to gravity 


I mass moment of inertia of one rotor blade 

about horizontal hinge 

K dimensionless coefficient (mB^R^2/2El) 

dL aerodynamic lift on blade element at radius r 


M bending moment In blade at radius r 

Mi» m 2» ^3 moments corresponding to the deflection 
J functions y^, yg, and y^ 

M- flexible-blade bending moment as defined In 

1 appendix B 


“r 


bending moment In blade at radius r (blade 
assumed to be a rigid body) 



weight moment of blade about horizontal 
hinge 
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m 

n 

P 

P 

Po 

Q 

R 

R Q 

r 

S 

t 

U T 

Up 

x 


mass of blade per unit length between 
hinge and radius BR 

an arbitrary integer 

any integer greater than zero and less than 
or equal to n 

tension in blade at radius r 

tension in blade at horizontal hinge 

mass of blade tip between radius BR and 
radius R 

blade radius 

distance from center of rotation to center 
of gravity of mass Q 

radius of blade element 

shear in the blade st radius r 

time 

velocity component at blade element perpen- 
dicular to blade span axis and parallel 
to rotor disk 

velocity component at blade element perpen- 
dicular both to blade span and to 

ratio of blade -element radius to BR 


y deflection of blade element at radius r, 

referred to rigid-blade position 


y 1? y2, yx deflection functions entering into the 

general equation for y 


blade flapping angle 

angle between plane perpendicular to axis of 
rotation and line connecting horizontal 
hinge with blade element at radius r 


P 

P* 
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^ slope of tangent to blade at radius r, 

referred to plane perpendicular to axis 
of rotation 

\|/ blade azimuth angle, measured In direction 

of rotation from down-wind position 

8 blade pitch angle 

0 angular velocity of rotor 

jx ratio of component of forward speed in plane 

perpendicular to axis of rotation to OR 

\ ratio of axial inflow velocity through rotor 

to OR 

y mass constant of rotor blade 

air density 
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Figure 1 Comparison of bending-moment analyses. 
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Figure 2.- Effect of tip-loss correction on bending moments 

computed by exact method . 
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Fig. 3 
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Figure 3.- Stresses due to combined centrifugal 
tension and bending in plane of flapping . 
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